The Imprint of Lithium Recombination on the 
Microwave Background Anisotropics 



Matias Zaldarriaga 1 and Abraham Loeb 2 
ABSTRACT 

o 

\ Following Loeb (2001), we explore the imprint of the resonant 6708A line 

opacity of neutral lithium on the temperature and polarization anisotropies of the 
cosmic microwave background (CMB) at observed wavelengths of 250-350//m (0.9-1.2 
THz). We show that if lithium recombines in the redshift range of z = 400-500 
■ as expected, then the standard CMB anisotropies would be significantly modified 

in this wavelength band. The modified polarization signal could be comparable to 
the expected polarization anisotropies of the far-infrared background on sub-degree 
^ angular scales (£ ^ 100). Detection of the predicted signal can be used to infer the 

■^j- ! primordial abundance of lithium, and to probe structure in the Universe at z ~ 500. 
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The latest measurements of the anisotropies in the cosmic microwave background (CMB; see 
Halverson et al. 2001, Lee et al. 2001, Netterfield et al. 2001) imply that the days of "precision 
cosmology" have already arrived^. Future ground and balloon based experiments, in combination 
with the satellite missions MAP[] in 2001 and PlanckP] in 2007, will test current theoretical models 
to a sub-percent precision at photon wavelengths ^ 500/iin. 

However, at far-infrared wavelengths of <J 350/im, Loeb (2001) has recently shown that 
neutral lithium can strongly modify the CMB anisotropy maps, through absorption and re-emission 
at its resonant 6708A transition from the ground state. Lithium is expected to recombine in the 
redshift interval z ~ 400-500 (Palla et al. 1995; Stancil et al. 1996, 1998). Despite the exceedingly 
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low lithium abundance^ produced in the Big Bang, the resonant optical depth (Sobolev 1946) 
after lithium recombination is substantial, 

r T , T [X(z)} = 1.00 1 , T (z) ( Xu , n \ ( ±±1) 32 . (i) 

v n Jhay ' V3.8 x 10- 10 J V 500 J v ' 

for an observed wavelength of \(z) = [6708Ax(l + z)] = 335.4//m x [(1 + z)/500]. Here, 
X Li ~ 3.8 x 10 -10 is the latest estimate of the lithium to hydrogen number density ratio (Buries 
et al. 2001), and f ul {z) is the neutral fraction of lithium as a function of redshift (Palla et al. 
1995; Stancil et al. 1996, 1998). Loeb (2001) argued that resonant scattering would suppress the 
original anisotropies by a factor of exp(— r LiI ), but will generate new anisotropies in the CMB 
temperature and polarization on sub-degree scales [l ^ 100), primarily through the Doppler 
effect. Observations at different far-infrared wavelengths could then probe different thin slices of 
the early universe. 

In this paper, we calculate in detail the effect of neutral lithium on both the polarization and 
temperature anisotropies of the CMB. §2 describes the modifications we have made to the standard 
code, CMBFAST (Seljak &; Zaldarriaga 1996), in order to calculate these anisotropies. In §3, we 
describe our results and compare them to the foreground noise introduced by the far-infrared 
emission from galaxies and quasars. Finally, §4 summarizes the main conclusions of this work. 

Throughout the paper we adopt the LCDM cosmological parameters of Oo ~ 0.25, Q\ = 0.7, 
Qb = 0.05, and Hq ~ 70 km s _1 Mpc -1 , and units of c = 1. 



2. Method of Calculation 

In order to compute the temperature and polarization fluctuations induced by lithium 
scattering, we complement the standard Thomson opacity (f T ) in CMBFAST by a new component 
(f LiI ) which is assumed to have a narrow Gaussian shape in conformal time, 

f = f T (7?)+f m (r),v) 

= an e a T + T LiI == , (2) 

where, r denotes optical depth, an overdot denotes a conformal time derivative, a is the expansion 
factor, rj is conformal time (a dr] = dt), n e is the number density of free electrons, and gt is 
the Thomson cross-section. The lithium opacity depends both on time and observed photon 
frequency, v. We characterize this opacity by three parameters: the total optical depth r LiI given 



6 Note that by the redshift of interest all the 7 Be produced during Big-Bang nucleosynthesis transforms to 7 Li 
through electron capture, since 7 Be starts to recombine well before 7 Li, owing to its significantly higher ionization 
potential. 
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by equation (|l|), and the central time and width (a ri ) of the Gaussian peak. The peak is 

set to a redshift (1 + z) = (A/6708 A) , where A = cjv is the observed wavelength. The actual 
(non-instrumental) width of the Gaussian is expected to be [Az/(1 + z)\ ~ 3 x 10~ 5 (Loeb 2001). 
Here we adopt a^jr) ~ 10~ 2 , as dictated by a fiducial detector's band-width. 

The lithium opacity introduces a frequency dependence to the CMB anisotropics. Next we 
show that different frequencies couple only through the net drag force they provide on the baryons. 
Our discussion follows the notation of Ma & Bertschinger (1995). 

We consider a density perturbation of comoving wave-vector k. The photon distribution 
function can then be expanded to first order in the perturbation amplitude as, 

oo 

*(k,n, !/,»;) = ^(-i)'(2/ + l)^(fc,^^(k-fi). (3) 

1=0 

where f (u) = 2/h 3 (e hu / kT ° - 1), and where / is a function of perturbation wave- vector (k), 
time (rj), and frequency (u) and propagation direction (n) of the photon. In the second line of 
equation (||), we have expanded the angular dependence of the distribution function in Legendre 
polynomials. 

The distribution function /(k, n, z/, rj) satisfies the Boltzmann equation. Since the Thomson 
cross-section is independent of frequency, the standard approach integrates the Boltzmann 
equation over photon frequencies and uses only one Boltzmann hierarchy to evolve the photon 
distribution function. In our case, the frequency dependence of the lithium cross-section implies 
that one should solve the hierarchy for different frequencies. Since the frequencies are coupled, one 
needs to follow all frequencies simultaneously, as done by Yu et al. (2001) for Rayleigh scattering. 
The coupling between different frequencies does not originate directly from the scattering term 
because lithium scattering does not change the photon frequency. As we show next, the coupling 
arises from the drag force on the baryons. Binary particle collisions allow the baryons to behave 
as a single fluid (Loeb 2001) which is subject to the sum of the forces applied by photons at all 
frequencies. In the limit of no drag force, the different frequencies decouple. 

Next, we define the relative density contrast (<5 7 ) and velocity divergence of the photon fluid 
(# 7 ) as, 

<5 7 = (o 4 p CMB ) -1 / d 3 v vf (v) •® (k,v,r]), 
3 

Se^u) = (a 4 ^^)" 1 -fc#i(/s,i/,77), 

e 1 = Jd 3 u vf (v) se^v). (4) 

With these definitions we can write the equation for the velocity divergence of the baryon fluid 
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where c s is the sound speed of the baryons, p CMB and p b are the mean energy densities of the 
CMB and the baryons, and 5 b is the baryon overdensity. The last term in the square brackets of 
equation (||) is responsible for the coupling between different photon frequencies. However, Figure 
1 of Loeb (2001) implies that we may ignore the drag force on the baryons (due to either Thomson 
or lithium scattering) at the redshifts where the lithium opacity becomes important (z <^ 500). 
Consequently, we may solve the Boltzmann equation separately for each photon frequency by 
explicitly neglecting the transfer of momentum from the photons to the baryons due to lithium 
scattering. Even though momentum conservation is not strictly satisfied in this approach, the 
remaining correction is expected to be negligible. 

The last change that we introduce to CMBFAST involves polarization. The cross-section for 
lithium scattering has a different dependence from Thomson scattering on both scattering angle 
and Stokes parameters. The scattering matrix for (Jii, I±, U) (where the parallel and perpendicular 
directions are defined relative to the scattering plane) can be decomposed into two parts, 



where is the scattering angle. The first term is the usual Thomson (or Rayleigh) scattering 
matrix multiplied by a factor E±, while the second term, which is proportional to E2, does not 
generate polarization and is isotropic in angle. The two amplitudes depend on the quantum 
numbers of the resonant states, and satisfy E\ + E2 = 1 (Chandrasekar 1960). For the transition 
between the ground state (2S) and first excited state (2P) of lithium, we get E\ = 1/3 (Hamilton 
1947; Chandrasekar 1960). Interference between the 2 2 S-2 2 P® /2 transition and the 2 2 S~-2 2 P° /2 
transition (Stenflo 1980) has a negligible effect on the polarization. This follows from the fact that 
t l .j ~ 1; and for a given lithium atom, a photon will likely scatter when its frequency is separated 
from the line center by less than the natural width (37 MHz/4-7r), which is much smaller than the 
frequency separation between these transitions (10 GHz). 



The significance of the new opacity component can be assessed from the visibility function 
T(r/). This function provides the probability distribution for the time of last scattering of the 
photons observed today at a conformal time 770, 
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Results 
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Fig. 1. — Visibility function for models with r LiI = 0.5 and LCDM. We show three examples, 
corresponding to resonant scattering at different redshifts, z = 400, z = 450 and z = 500. The 
present-day value of the conformal time is r]o = 1.39 x 10 4 Mpc. 

r(rj) = [ V ° drfrti). (7) 

In Figure [l] we show the visibility functions for some of the models we consider later. 

The observed anisotropics have two separate contributions, one from the standard last 
scattering surface at hydrogen recombination (decoupling), which is suppressed by e _T Ln, and a 
second new contribution that is generated by lithium scattering at lower redshifts. In the following 
sub-sections we will characterize this new contribution to the temperature and polarization 
anisotropies. 



3.1. Temperature and Polarization Anisotropies 

Figures [2] and ^ show the predicted power spectra for the temperature and polarization 
anisotropies of the CMB at an observed wavelength of A = 335.4/xm, corresponding to lithium 
scattering at z = 500. The Stokes parameters are measured in /j>K. The figures compare the 
spectrum of fluctuations in standard LCDM (no lithium scattering) with two other models, each 
having a peak in f LiI at a redshift of z = 500, but with a total optical depth of either r LiI = 0.5 or 
t ui = 2.0. Since at long wavelengths, A ^ 500/um, the LCDM fluctuations are not altered, precise 



Fig. 2. — Temperature power spectra for the standard calculation (LCDM model), and two other 
models with lithium optical depths of r LiI = 0.5 or r LiI = 2.0 at z = 500, for observations in a 
narrow-band around 335/im. The width of the f LiI Gaussian in Eq. (g) is cr^/rj = 0.01. 



mapping of these fluctuations by the MAP or Planck satellites will provide a reference power 
spectrum against which the lithium distortion can be measured. 

Figure ^ shows that for r LiI = 0.5 the small scale fluctuations are dominated by the suppressed 
anisotropies from recombination, resulting in a power spectrum which is similar in shape to that 
of the primary anisotropies but suppressed in amplitude. However, the r LiI = 2 case is very 
different. Here, the anisotropies are actually larger for many Us than those expected without the 
lithium scattering, having a different functional dependence on I than the standard case. The 
e -T Lii suppression of the original anisotropies is sufficient to make them sub-dominant relative to 
the newly generated anisotropies at z = 500. 

There is an interesting difference between the primary anisotropies and those created by 
lithium scattering. In order to explain it, we introduce the integral solution for the temperature 
anisotropies, 



AT 



no 



+ (k, n, r, ) = / dr, (5 7 /4 + V + A»&)Tfa) + (</> + ^)e~ T ^ 



(8) 



where [i is the cosine of the angle between the wave-vector and the direction of observation 
and (j) and V are the two gravitational potentials defined by the perturbed metric, 
ds 2 = a 2 (rj)[— (1 + 2ip)dr]' 2 + (1 — 2<j))dxidx 1 ]. While the contribution from recombination 
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is dominated by the monopole term, (5 7 /4 + if)), the lithium anisotropies are dominated by the 
peculiar velocity term (oc v^) on most scales. We illustrate this in Figure |3| where we show the 
monopole and velocity contributions to the anisotropies in both the standard LCDM model and in 
a model that has a large optical depth r LiI = 10 at z = 500. We chose such a large optical depth 
in order to suppress the original contribution from decoupling. The figure clearly shows that the 
anisotropies are dominated by the monopole term for the standard LCDM model for almost all 
values of I, while the opposite is true for the r LiI = 10 model. The new anisotropies are dominated 
by the monopole term only at very low multipoles, I <J 20. 

We can easily explain why the monopole no longer dominates for the lithium contribution. 
After recombination the monopole term decays by the free streaming of the photons while the 
velocity of the baryons continues to grow as they fall into the dark matter potential wells. For a 
perturbation mode of wave- vector k, the monopole term at conformal time r] after recombination, 
f] > Vrec, is approximately given by, 

(&y/4 + VOfa) = (V 4 + ^)farec)jO [HV ~ Vrec)] , (9) 

where jo(x) is the spherical Bessel function. Equation @ shows that the monopole term is small 
for k(rj — 7] TCC ) 3> 1 because of the decay in the Bessel function when its argument is large. Figure 
[l] shows that the new peak of the visibility function occurs at r] ~ 500 Mpc while r/ TCC ~ 300 Mpc. 
We can translate the spatial wavenumber k to angular scale using the conformal distance to the 
new peak in the visibility function, d = ~~ V) ~ Vo- We find that k(rj — r/ TCC ) ~ 1.6 x 10 — 2 Z, 
which explains why the monopole term is suppressed for I 3> 60. 

While the monopole term decays between recombination and the lithium scattering surface, 
the velocity grows and thus produces anisotropies that are larger than those generated at 
decoupling in the r LiI = 2 case. 

The physics of the polarization anisotropies is different from that of the temperature 
anisotropies. Since polarization is generated by the quadrupole moment, there are two competing 
effects that need to be considered. On the one hand, the quadrupole anisotropies are small at 
recombination, since they are suppressed relative to the velocity fluctuations by a factor k drj, 
where 5r] is the width of the last scattering surface at recombination (Zaldarriaga & Harari 1995). 
In the new scenario the quadrupole is able to grow during the free streaming period between 
recombination and z ~ 500. This naturally leads to an increase in the polarization signal. The 
same effect increases the polarization anisotropies on large scales in models with a substantial 
optical depth to Thomson scattering after the universe reionizes (Zaldarriaga 1997). On the 
other hand, due to the nature of resonant -line scattering (Hamilton 1947; Chandrasekhar 1960), 
only 1/3 of the cross-section generates polarization out of this quadrupole, while 2/3 produces 
unpolarized radiation (E\ = 1/3, E2 = 2/3). Although the quadrupole is bigger at z ~ 500 
than at recombination, the newly generated polarization is not as large. For example, Figure || 
shows that even for r UI = 2, the polarization at high multipoles, I ~ 1000, is dominated by the 
suppressed signal from decoupling. 
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Fig. 3. — Velocity and temperature contribution to anisotropies for LCDM and for a model with 
r Lii = 10 at z = 500. The curves that approach a finite value at low I are the temperature 
contributions, and the curves that approach zero at low I describe the velocity contributions. 

3.2. Cross correlations and Difference Maps 

Let us next consider two temperature anisotropy maps, obtained from observations at different 
wavelengths, T a {6) and T b (6). We define the following spectra in multipole space, 

{t?t$) = s w cr, 

(T^) = 8 lv Cf, 

(T?T$) = 5 lv Cf, (10) 

with being the Kronecker delta. Similar expressions can be written for E and -B-type 
polarization. We then define the cross-correlation coefficient, 

CCi = . 1 . (11) 

If CC\ — > 1 across a range of l's, then the two maps are scaled versions of each other over that 
range. 

Figures [5] and |6| show the correlation coefficients between maps of the anisotropies at different 
observed wavelengths, labelled by the redshift of the peak in the visibility function. They 
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Fig. 4. — Polarization power spectra. Models are the same as in Fig. |2[ 

correspond to the correlation coefficient in three cases: (i) between the primary anisotropics and 
those measured at A = 335.4 /im (z = 500); (ii) between measurements at A = 335.4 /im and 
A = 301.8 fjan = 450); and (iii) between A = 335.4 ,um and A = 268.3 //m (z = 400). 

The temperature correlation coefficient between the original map of the primary anisotropies 
and a map at 335.4 ^m can be understood as follows. For low I, CCi « 1 because the anisotropies 
are produced by long wavelength (small k) modes and the difference in conformal time between 
recombination and z ~ 500 (which we denote by d rec -5oo) is smaller than the perturbation 
wavelength , i.e. (k d rec _5oo << 1). The fact that there are two different scattering surfaces does 
not make a difference for these modes. 

On small scales, the situation is more complicated. The component of the anisotropies from 
decoupling, which is suppressed by a factor e~ r Ln , drives the cross-correlation coefficient to unity, 
because its suppression results merely in rescaling. However, the newly generated anisotropies are 
uncorrelated with those coming from decoupling if they are produced by wavelengths that are 
smaller than d rec _5oo, i-e. if k d rec -500 ~ (^rec-soo/^rec)^ ~ 1-6 x 10 2 Z ^> 1. The newly generated 
anisotropies tend to drive the cross-correlation coefficient to zero on small scales. 

In Figure [| we also show the cross-correlations between maps at observed wavelengths for 
which f LiI peaks at z = 400, z = 450, and z = 500. For a given angular scale, the anisotropies 
generated by the second peak of the visibility function are different only if the conformal distance 
between the two peaks is much larger than the wavelength of the perturbation producing the 
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Fig. 5. — Correlation coefficients for temperature anisotropies at three different wavelengths 
corresponding to lithium resonances at z = 400 (268/um), z = 450 (302/mi) and z = 500 
(335//m). The dotted lines show the correlation coefficient between the standard (long wavelength) 
anisotropies and those with lithium scattering at z = 500; the solid line shows the correlation 
between the anisotropies for 302 /im and 335/im (z = 500 and z = 450); and the dashed lines 
correspond to 335/mi and 268/mi (z = 500 and z = 400). The upper panel shows the case of 
r LiI = 0.5 and the lower panel shows the case of r LiI = 2.0. 

anisotropies on that scale. This explains why the departure from CC\ = 1 occurs on smaller scales 
for the correlation between z = 500 and z = 450 than for that between z = 500 and z = 400. 
Also, if the anisotropies in the two maps have a substantial contribution from decoupling, then 
the correlation coefficient will not approach zero even on small scales where the newly generated 
anisotropies in the two maps are uncorrelated. The significance of this primary contribution from 
decoupling is suppressed when r LiI increases, as illustrated by the r LiI = 0.5 and r LiI = 2.0 panels 
of Figure ||. 

The correlation coefficient for polarization which is shown in Figure |6|, can be explained along 
similar lines. On small scales, the polarization anisotropies are dominated by the suppressed 
contribution from recombination, and so CCi — > 1. The dominance of the new anisotropies over 
the primary ones at I ~ 200, makes the correlation coefficient between primary and z = 500 deviate 
away from unity around that scale. However, among neighboring frequencies the cross-correlation 
approaches unity near I ~ 200, because the distance between the peaks of the visibility function is 
not sufficient to decorrelate the contributions from the relevant k modes. 
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Fig. 6. — Correlation coefficient for polarization. The upper panel shows the case of r LiI = 0.5, and 
the lower panel shows the case of r LiI = 2.0 

In Figures and || we show the difference power spectra for two maps at different observed 
wavelengths. We define <5T(0) = T a {6) - T b {6) and 

{5Tm,) = 5 lv Cf m . (12) 

A similar expression can be written for the polarization. The difference spectrum provides the 
new signature due to lithium, since the MAP or Planck satellites will measure with high precision 
the anisotropies at long photon wavelengths. 

The power spectra of the difference maps, shown in Figure ^ for the temperature and Figure 
H for the .E-type polarization, are consistent with our interpretation of the cross-correlation 
spectra between maps. On large scales, the difference between the temperature maps becomes 
very small. The power in the Z500-Z400 difference map peaks at I ~ 450, while the power for the 
recombination- ZbOO difference peaks at I ~ 250. The amplitude of the difference increases as the 
optical depth increases. When r ha is large, the difference spectrum for the recombination- Z 500 
polarization almost coincides with that of the primary anisotropies, as a result of the fact that the 
primary anisotropies dominate over the newly generated ones. 
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Fig. 7. — Difference power spectra for temperature anisotropies. In the upper panel r LiI = 0.5, and 
in the lower panel r LiI = 2.0. 

3.3. Comparison to the Far Infrared Foreground 

The main difficulty in measuring the lithium imprint on the CMB anisotropies is the 
contamination by the far-infrared background (FIB). There has been no detection to date of the 
anisotropies in this background, and so we have to rely on theoretical estimates (Haiman & Knox 
2000; Knox et al. 2001). 

In order to properly combine the contributions from the CMB and FIB anisotropies we 
express intensities in terms of the equivalent Rayleigh-Jeans temperature, T R j (in fiK). We start 
by comparing the temperature anisotropies. The total fluctuation amplitude is given by, 

AT A rpCMB , at^FIB a rp rpCMB (ATrj\ ^FIB /ATrjN 
Ai RJ =Ai RJ +AJ RJ =AJ RJ = i RJ — +i RJ — . (13) 

V J RJ / CMB V 1 RJ /FIB 

The ratio between the CMB intensity and the central value for the inferred intensity of the 
FIB (Fixsen et al. 1998; see the central dot-dashed curve in Figure 2 of Haiman & Knox 2000) is 
of order unity for a lithium scattering redshift z ~ 500, 

t rt MB ( 500 \ f m ( 500 \1 
As noted by Loeb (2001), the temperature fluctuations in the Wien tail translate to intensity 
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Fig. 8. — Difference spectra for polarization. In the upper panel r LiI = 0.5 and in the lower panel 
r m = 2.0. 

fluctuations, Deltal u (in erg s~ 1 cm~ 2 sr _1 Hz~ 1 ), of a much larger contrast, 

fAT KA _ AI U fdlnI A AT ( AT ^ x /_500_N AT 





AT 
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\dhxTj 


~T~ ~ 




-T 



V Irj /cmb ^ WlnT/ T V&T/ T \1 + z J T 

where we have substituted I V (T) oc exp{—hv/kT) and T = 2.725 K (Mather et al. 
1999). The anisotropy amplitude shown in Figure ^ depends on r LiI , but roughly implies 
(AT R j/Trj)cmb ~ 3 x 10" 4 [500/(1 + z)\. Haiman & Knox (2000) and Knox et al. (2001) estimate, 
(ATrj/Tjjj)fib = 0.05-0.1. The FIB anisotropics peak at an I of a few hundred but the peak is 
very broad. The anisotropies in the FIB are relatively large since they originate from clustering of 
sources at low redshifts, z ~ 1. 

We conclude that if 50% of the lithium ions recombine at z ~ 500 then 

a^I' b „ 3 x 10 -s ( my f 15J8 A _ 500 N 1 _ 



at|^ b vi + ^y ' i v 1 + 2! 

Since the CMB contribution is sub-dominant, it is essential to exploit the different frequency 
dependence of the FIB and CMB anisotropies in order to subtract the FIB contribution with high 
precision. This might be possible on small angular scales, where the temperature anisotropies 
generated by lithium scattering at different redshifts are uncorrelated, as indicated by Figures |5| 
and f?]. Also, since the FIB is produced by point sources, observations with high angular resolution 
can resolve the sources and remove them individually. 
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Next we consider polarization. While the intensity anisotropics of the FIB are dominated by 
the clustering of the sources rather than by Poisson fluctuations, the situation is different for the 
polarization anisotropies. We assume that each source is polarized with a degree of polarization 
e, and that the polarizations of different sources are uncorrelated. Under these assumptions, it is 
easy to show that both the E and _B-polarization spectra share the same amplitude of 

e 2 /2 x C|r s ° n , (17) 

where C^° lsson is the Poisson contribution to the temperature anisotropies. The Poisson 
fluctuations can be calculated based on the SCUBA source counts, and are shown as the light 
solid line in Figure 3 of Haiman & Knox (2000). 

The Poisson component of the FIB anisotropies is approximately a factor of ~ 100 smaller in 
power than the total temperature anisotropies at I ~ 100-200. If we adopt an average value of 
e = 0.4% (Jones 1993), then the ^-polarization is lower by a factor (0.4%/\/2) x 0.1 = 2.8 x 10" 4 
than the intensity fluctuations. 

On the other hand, Figures |2| and |] imply that the CMB is approximately 10 % polarized for 
I ;> 200. We therefore get that as long as lithium recombines around z ~ 500, 

AES¥ b ( e \ ( 500 \ 2 f / 500 \1 

3HTr»l-lX ^7 (ttT «xp 15.78 1 --- . (18) 



AE™ \0A%J\l + zJ *\ V 1 

A property which could be very useful in experimental attempts to isolate the signal is 
that the polarization signature due to lithium has only an E-type component, while the FIB 
polarization has an equal amplitude in both the E and B— type components. The 5— type 
polarization could therefore be used to monitor the FIB contamination and could play an essential 
role in the subtraction of the FIB. 

Depending on the nature of the sources responsible for the FIB and their luminosity function 
it may eventually become possible to resolve most of the FIB through high-resolution observations 
at different wavelengths. This approach is used, for example, in observational studies of the 
Sunyaev-Zeldovich effect, where much of the contribution from discrete foreground sources is 
subtracted out through deep, high-resolution observations at either radio or optical-infrared 
wavelengths. At the present time there are no available source counts in the wavelength range we 
consider here. Closest in wavelength are source counts from the SCUBA instrument (see, e.g. Fig. 
2 in Borys et al. 2000). If most of the FIB could be resolved, the task of detecting the effect of 
lithium would become easier as the overall level of contamination would be drastically reduced. 
Future studies of the FIB will determine whether this reduction is feasible. 
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4. Conclusions 

We have shown that if more than half of the lithium ions recombine by z ~ 500, then the 
temperature and polarization anisotropies of the CMB would be strongly altered at an observed 
wavelength of 335^m (see Figs. ^ and For high multipoles I 3> 10, the change is dominated by 
two contributions: (i) the Doppler anisotropies induced at the sharp lithium scattering surface; 
and (ii) the uniform exp{— r LiI } suppression of the primary anisotropies which were generated 
at hydrogen recombination (decoupling). Maps taken at wavelengths that are different by only 
~ 10% are expected to have significant differences (see Figs. ^ and|8|). 

The above signals are superimposed on top of the far infrared background (FIB). Our 
estimates imply that the lithium imprint on the CMB polarization should be comparable to that 
provided by the FIB (Eq. p^). Detection is more difficult for the temperature anisotropies (Eq. 

HP- 

The wavelength range we explored overlaps with the highest frequency channel of the Planck 
mission (352/iin) and with the proposed balloon-borne Explorer of Diffuse Galactic Emissions]] 
(EDGE) which will survey 1% of the sky in 10 wavelength bands between 230-2000/im with a 
resolution ranging from 6' to 14' (see Table 1 in Knox et al. 2001). In order to optimize the 
detection of the lithium signature on the CMB anisotropies, a new instrument design is required, 
with multiple narrow bands (AA/A <^ 0.1) at various wavelengths in the range A = 250-350/iin. 
The experiment should cover a sufficiently large area of the sky so as to determine reliably the 
statistics of fluctuations on degree scales. In order to minimize contamination from the FIB, the 
detector should be sensitive to polarization. For reference, the experiment should also measure 
the anisotropies at shorter wavelengths where the FIB dominates. In order to detect the effect of 
lithium, high signal-to-noise maps of the primordial CMB should be made for the same region of 
the sky. Most likely, those maps will become available from future CMB missions such as Planck. 
A strategy for eliminating the contribution from the brightest FIB sources may also be needed. 

The resonant optical depth depends sensitively on the primordial lithium abundance and the 
recombination history of lithium. More detailed calculations of lithium recombination will be done 
in a forthcoming paper (Dalgarno, Loeb, & Stancil 2001). Detection of the lithium signature will 
also allow to calibrate the primordial lithium abundance, which is a sensitive indicator of the mean 
value and the dumpiness in the baryon abundance during Big Bang nucleosythesis. The lithium 
abundance in nearby stars is subject to large astrophysical uncertainties (Buries et al. 2001, and 
references therein). We note that values of the lithium opacity which are higher than the ones we 
have used, are potentially possible. As an extreme example, lithium abundance values as high as 
Xlii ~ 10 -8 were suggested by models of inhomogeneous Big-Bang nucleosynthesis (Applegate & 
Hogan 1985; Sale & Mathews 1986; Mathews et al. 1990), 



7 http: / /topweb. gsfc.nasa.gov 
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The lithium signature on the CMB anisotropies provides the only direct probe of structure 
in the universe at a redshift z ~ 400-500. This redshift marks the beginning of the "dark ages" 
which end only after the first generation of galaxies form at z ~ 20 (see review by Barkana & Loeb 
2001). 
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